Fermi surface nesting and charge-density wave formation in rare-earth tritellurides 
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The Fermi surface of rare-earth tri-tellurides (7?Te3) is investigated in terms of the nesting driven 
charge-density wave formation using positron annihilation and first-principles LMTO calculations. 
Fermi surface nesting is revealed as a strong candidate for driving charge-density wave formation in 
these compounds. The nesting vector obtained from positron annihilation experiments on GdTe3 
is determined to be q = (0.28 ± 0.02, 0, 0) a*, (a* = 27r/a), in excellent agreement with previous 
experimental and theoretical studies. 

PACS numbers: 71.45.Lr,71.18.+y,78.70.Bj 



I. INTRODUCTION 

The existence of charge-density wave (CDW) ground 
states in low dimensional metals has attracted a great 
deal of interest over many years (see for example 
It is well known that instabilities in the Fermi surface 
(FS) due to particular features of the electronic band 
structure (e.g. nesting features, van Hove saddle points) 
can lead to the emergence of new ground states such as 
spin-density Q or charge-density waves Examples 
of systems where this is thought to be the case include 
the transition metal dichalcogenides [[J, pure Cr Q and 
organic charge-transfer salts |2j|. While the underlying 
electronic band structure of a system whose new period- 
icity is commensurate with the lattice is well understood, 
the behavior in cases where the periodicities are incom- 
mensurate is much less clear 

Motivated by recent angle-resolved photoemission 
(ARPES) studies aimed at clarifying this issue we 
have combined first-principles electronic structure calcu- 
lations and positron annihilation experiments to establish 
the FS topologies of a model CDW system, namely the 
rare earth tritellurides. 

In systems whose FS topologies include large parallel 
faces spanned by a (nesting) vector q, there will be an in- 
stability in that system towards opening up gaps on these 
faces by introducing some new periodicity and hence new 
superzone boundaries. A CDW is an example of the re- 
sults of such an instability. Although in general the CDW 
will not be commensurate with the lattice (since its pe- 
riodicity is dictated by the nesting vector), many CDWs 
experience "lock-in" to a suitably commensurate period, 
thereby lowering the strain energies associated with the 



distortion. 

The rare-earth tritelluride series, i?Te3 (R = Y, La - 
Tm), crystallises in the layered NdTe3 structure and 
consists of sheets of inert i?Te slabs sandwiched between 
a double-layer of square planar Te sheets (see Fig. QJ. 
It should be noted that the crystollagraphic &-axis is the 
long axis. Resistivity measurements have revealed a large 
anisotropy (p a /pb) across the series, as much as ~ 5000 in 
the case of SmTe3 01 • The electron diffraction results of 
DiMasi et al. showed strong satellite peaks, which were 
interpreted as stemming from the periodic lattice dis- 
tortions associated with the presence of an incommen- 
surate CDW They report a range of wavevectors 
(q « 0.27-0.31 a*, a* = 2n/a) (i.e. around (2/7) a*) for 
i?Te3, where the CDW was found to be stable for all the 
rare-earths, R (from La to Tm), and occurred for most 
materials along both the a and c axes. In some samples, 
CDW formation was only observed along one of these 
directions, although this was found to be highly sam- 
ple dependent. ARPES studies 0, have successfully 
mapped the gap anisotropy, finding the FS to be very 
significantly gapped (~ 200 meV) along r — Y (T — Z in 
their coordinates); this would correspond to a tempera- 
ture well beyond the samples' melting points. Further, 
the ARPES studies have confirmed the magnitude of the 
CDW wave- vector to be q « (2/7) a*, equal to that ob- 
served by DiMasi et al. 0|. The RTe^ system thus 
presents an unprecedented opportunity to study the cor- 
relation between the FS topology (through its nesting) 
and the CDW period. 

The rare-earth ditellurides RTe2 (Fig. are a re- 
lated material, crystallising in the layered tetragonal 
anti-Cu2Sb structure 0|. In these materials, electron 
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diffraction reveals satellite peaks corresponding to a com- 
mensurate superstructure characterised by a (2 x 1 x 1) 
supercell, attributed to CDW formation Calcula- 
tions of the electronic structure of this system will be 
presented alongside those for RTe^ in order to highlight 
the similarities and differences. 

Here we present the FS of GdTe3 as mea- 
sured by 2-dimensional angular correlation of electron- 
positron annihilation radiation (2D-ACAR) in conjunc- 
tion with first-principles electronic structure calculations 
of lutetium di- and tri-tellurides. 



II. ELECTRONIC STRUCTURE 
CALCULATIONS 

The electronic band structure was calculated for both 
lutetium di- and tritelluridc using the linearized muffin- 
tin orbital (LMTO) method within the atomic sphere ap- 
proximation (ASA) , including combined-correction terms 
151. Lu was chosen to avoid the complications associated 
with the description of /-electrons within the local den- 
sity approximation. The parameters of the crystalline 
structure were 4.55A (b/a = 2.037) and 4.34A (b/a = 
5.93) respectively, corresponding to those found for 
LaTe 2 in Ref. Q and SmTe 3 in Ref. The results 

of the calculation were found to be relatively insensitive 
to changes in the lattice parameter of ~ 5%, allowing us 
to address all the different rare-earth systems 01 • All 
calculations used a basis of s, p, d and / states, and 
self-consistency was achieved at 1568 and 1176 k-points 
respectively in the irreducible Brillouin zone (BZ). 

Although the tritcllurides are found to have very 
slightly orthorhombic unit cells (a / c), a square-based 
cell (a — c) is adopted (since the effects were found to 
be insignificant). It is worth highlighting that there will 
still be in-plane anisotropy in the tritcllurides even for 
a = c, caused by the RTe layer orientation, which will be 
reflected in the electronic structure. 



A. Lutetium ditelluride (LuTe2) 

The electronic band structure and FS for LuTe2 is 
shown in Figs.|2and|21respectively. As expected, there is 
little electronic interaction between layers in these com- 
pounds, expressed in the small dispersion along [010] (e.g. 
r — Y and N — R). We observe three bands to cross 
the Fermi level (-Ep), principally due to 5p orbitals of 
the planar Te layer. Two of these bands result in the 
diamond-shaped FS, also predicted by previous calcu- 
lations [Til \wi Hsj . These are the sheets which show 
the nesting. Two further FS sheets are observed in the 
present calculation, forming small hole pockets centered 
at Y (see Fig. HJ). As may be seen in Fig. |21 the bands 
responsible for these two sheets only just cross -Ep, and 
thus the presence or absence of these sheets is likely to 
be very sensitive to the details of the calculation (such as 



the description of the Lu 4/ electrons in the local den- 
sity approximation) . It is emphasised that these features 
are irrelevant to the nesting vector previously proposed 
by DiMasi et al. 0], corresponding to a commensurate 
nesting vector of q = (1/2) a*, shown in Fig. |3| our 
results are consistent with the existence of such a nesting 
vector. 



B. Lutetium tritelluride (LuTes) 

The LuTe3 band structure (Fig. HJ differs most no- 
tably from LuTe2 in the appearance of a bilayer splitting 
induced as a result of the coupling between the two Te 
planes. This is observed as a small splitting of the sheets 
which compose the FS (Fig. EJ) , and which has been re- 
ported ex peri mentally in ARPES studies for CeTe3 and 
SmTe 3 0,111. 

Fig. shows the nesting vector previously proposed 
for these materials, q ~ (2/7) a*, by DiMasi et al. [Til ], 
observed directly in ARPES studies E3 • It is clear 
from this figure that our band structure results strongly 
support such potential for nesting, further emphasised 
by the calculation of x(<?)j f° r Q along [100] (see Fig.EJ. 
This displays a clear, broad peak centered at q « 0.26 a* 
(the width of which arises due to the splitting of the FS), 
corresponding to nesting between the inner and outer 
diamond-shaped features (as shown in Fig. • 

Since the FS is composed of states associated with the 
Te layer, no appreciable difference is likely to be observed 
in the nesting features of -RTC3 for different rare earth 
atoms, indicated by the stability of the wavevector for 
different R observed by DiMasi et al. ^H, and so com- 
parisons may be made between the lutetium compounds 
used in the calculation and the gadolinium compounds 
used in the experiment (and, for that matter, the rare- 
earths in other studies). 



III. POSITRON ANNIHILATION 
MEASUREMENTS OF THE FERMI SURFACE 

The occupied momentum states, and hence the FS, can 
be accessed via the momentum distribution using the 2D- 
ACAR technique 0]. This well-established technique 
has recently been used by some of the present authors to 
determine the FS topology and identify nesting features 
in a wide range of systems 20] . A 2D-ACAR measure- 
ment yields a 2D projection (integration over one dimen- 
sion) of the underlying momentum distribution. Given 
that the RTe^, series is strongly two-dimensional, the dis- 
persion along the T — Y direction is small and thus the FS 
is strongly two dimensional; in this case one single pro- 
jection along [010] is sufficient to observe the topology 
of the FS. Application of the Lock-Crisp- West (LCW) 
prescription |2l| . where the p-space momentum density 
is 'folded' back into the first BZ, yields a translation- 
ally invariant k-space distribution which indicates the 
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level of occupancy of the states across the (projected) 
BZ. The FS topology expresses itself through changes in 
this occupancy according to whether a particular k-state 
is occupied or unoccupied for each band. The virtue of 
the 2D-ACAR technique in such studies is that it reveals 
directly the shape of the FS through the shape of the 
occupied and unoccupied regions. The effect of a charge- 
density wave gap on the momentum distribution will be 
to slightly smear the step associated with the presence 
of the FS; however, given that the intrinsic momentum 
resolution of the spectrometer is ~ 15% of the BZ, this 
additional smearing is negligible and hence the 2D-ACAR 
spectrum is insensitive to the presence of the CDW. The 
positron annihilation technique thus reveals the FS un- 
reconstructed by the CDW (i.e. the ungapped FS). 

The positron annihilation measurements were per- 
formed using the Bristol spectrometer. The single crystal 
samples were grown by slow cooling of a binary melt con- 
taining an excess of Te, as described previously by Ref. 
0- The GdTe3 sample was oriented to obtain a projec- 
tion along the [010] crystallographic direction, and ~ 150 
million counts were accumulated at a sample temperature 
of ~ 80K. A maximum entropy based deconvolution rou- 
tine was used to reduce the smearing caused by the finite 
instrument resolution [22j . 

The k-space occupation density, projected down the 
[010] direction, is presented in Fig. [7| alongside the oc- 
cupancy obtained from the LMTO calculation. Fur- 
ther measurements of TbTe3 were also performed at a 
lower statistical precision, revealing identical features to 
GdTe 3 , and will not be discussed any further. Excellent 
agreement is observed between the experimental and the- 
oretical BZ occupancy; the X electron pockets are clearly 
visible in the data, and agree in both size and shape 
with the theoretically determined pockets. Whereas in 
the ARPES studies of Ref. 12] the FS was observed to 
be gapped along T — Y (T — Z in their coordinates), in 
this study no such gapping is observed, the effect of the 
CDW being negligible when compared with the momen- 
tum resolution of the apparatus (as described above). 

As the LMTO calculation highlights, the Te bands 
that cross the Fermi level are degenerate for the ditel- 
lurides. However, owing to coupling between the two 
non-equivalent square planar Te sheets in the tritel- 
lurides, these split. Whilst this splitting has been ob- 
served (for example, 0, 0), it should be pointed out 
that our momentum resolution precludes observation of 
this small splitting of the FS in our experiment. Note, 
however, that the splitting, being smaller than the exper- 
imental resolution, is likely to enhance the observability 
of the nested Fermi surface. 



The interpretation of the data in terms of nesting 
provides compelling evidence in support of the role of 
the FS in driving CDW formation. Strong nesting fea- 
tures, identical to those outlined in the calculated FS, 
are clearly visible in the data, connecting large paral- 
lel surfaces of FS (Fig. 0), and revealing a nesting vec- 
tor (present in both [100] and [001] directions) connect- 
ing optimally nested portions of the experimental FS of 
q = (0.28 ± 0.02,0,0) a*. (The Fermi surface breaks 
have been determined from the experimental data by a 
"zero crossing" method [2^| , where the "raw" k-space oc- 
cupancy map has been subtracted from the "Maximum 
Entropy deconvoluted" k-space occupancy of the raw 
data. The zero-level contour reveals the breaks due to the 
diamond-shaped sheets, from which the nesting vector 
has been extracted.) This compares favorably with the 
superlattice reflections observed in electron diffraction 
experiments and with nesting vectors of q m (2/7) a* 
revealed via ARPES studies fHn%. 



IV. CONCLUSION 

In conclusion, we have presented the experimental 
ungapped Fermi surface of the high-temperature CDW 
compound GdTe3. Complementary band structure cal- 
culations on LuTe3 and the related compound, LuTe2, 
are also presented. Both the experimental and theoreti- 
cal FS topologies display excellent propensity for nesting, 
supporting the role of the Fermi surface in driving the 
formation of the CDW in these compounds. The nest- 
ing features are found to be in excellent agreement with 
previously proposed experimental an d theoret- 

ical [ll], HH US studies on the i?Te3 series, revealing an 
experimental nesting vector of q = (0.28±0.02, 0, 0) a*. 
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FIG. 1: (color online) The structure of (a) 7?Te2 and (b) 
RTe3. Note the two Te planes sandwiched between itTe layers 
in 7?Te 3 . 
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FIG. 2: The electronic band structure of LuTe2. 
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FIG. 3: (color online) The calculated Fermi surface of LuTe2 
(a) on a (010) plane through V with the proposed nesting 
vector illustrated (for clarity the Y pockets have been omitted 
here), (b) the FS sheets responsible for the nesting, and (c) 
the F-centered hole pockets. High symmetry points have been 
labeled in (c). The F-point is located at the center of the zone. 
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FIG. 4: The electronic band structure of LuTe3. 
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FIG. 5: (color online) The calculated Fermi surface of LuTe3 
(a) on a (010) plane through V with the proposed nesting 
illustrated, (b) the two diamond shaped sheets (split by the 
coupling between the two Te planes) which exhibit strong 
nesting, and (c) the two X pockets. The symmetry points 
are the same as those labelled in Fig. 
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FIG. 6: The generalized susceptibility x(o) calculated using 
the LMTO method for LuTe3 along the [100] direction, show- 
ing the broad peak around q ~ 0.26 a* due to nesting between 
the two diamond shaped sheets. 



9 




FIG. 7: (color online) (a) Projected occupation density ob- 
tained from the LMTO calculation for LuTe3, and (b) occu- 
pation density within the first BZ for GdTe3 measured by the 
2D-ACAR experiment. In both cases, lighter colors (shades) 
represent higher occupancy, and high symmetry points have 
been labelled in projection. 



